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Snapshots  of  High-Latitude  Electrodynamics 
Using  Viking  and  DMSP  F7  Observations 


G.  T.  Marklund,*  L.  G.  Blomberg,'  K.  Stasiewicz,*’®  J.  S.  Murphrbe,^  R.  Pottblbtte,® 
L.  J.  Zanbtti,®  T.  a.  Potemra  ®  D.  A.  Hardy/  and  F.  J.  Rich^ 


Simultaneous  observations  by  the  Viking  and  the  DMSP  F7  satellites  have  been  used  as  input 
to  a  new  method  to  obtain  snapshot  pictures  of  the  auroral  electrodynamics.  In  particular,  an 
"instantaneous”  global  eqiupotential  (or  convection)  pattern  is  calculated  from  distributions  of 
field-aligned  current  and  conductivity  which  are  qualitatively  consistent  with  the  Viking  atiroral 
imager  data  and  quantitatively  consistent  with  magnetic  field  and  particle  data  from  the  two 
satellites.  This  convection  pattern,  which  is  of  the  normal  two-cell  type,  with  a  weak  dusk  cell  and 
a  strong,  elongated  crescent-shaped  dawn  cell  (consistent  with  positive  interplanetary  magnetic 
field  By),  agrees  well  with  the  Viking  electric  field  data.  The  model  and  the  observed  potential 
profiles  agree  nicely  along  the  entire  Viking  orbit  except  for  two  intervals  above  acceleration 
regions  where  deviatioiu  are  to  be  expected  (due  to  parallel  electric  fields).  These  regions  are 
characterised  by  D -shaped  potential  minima,  upward  field-aligned  currents,  upgoing  ion  beams, 
and  relatively  intense  auroral  kilometric  radiation.  Thus,  the  model  results  are  consistent  with 
the  Viking  observations  not  only  on  a  global  scale  but  also  on  the  scale  of  the  auroral  acceleration 
regions.  The  corresponding  convection  in  the  magnetosphere  is  obtained  from  a  simple  projection 
to  the  equatorial  plane  of  the  deduced  two-cell  convection  pattern.  FVom  this  the  location  of  the 
plasmapause  is  inferred 


1.  Introduction 

In  our  efforts  to  better  understand  the  magnetosphere- 
ionosphere  system  a  major  problem  is  to  provide  an  accept¬ 
able  coverage  by  simultaneous  observatioas  within  the  large 
volume  to  be  studied.  The  high-latitude  ionosphere  acts 
in  many  respects  as  a  projection  screen  for  processes  tak¬ 
ing  place  in  the  outer  magnetosphere  and  boundary  regions. 
Much  can  therefore  be  learnt  from  studies  of  the  electro¬ 
dynamics  of  the  high-latitude  ionosphere,  which  is  a  rela¬ 
tively  small  and  rather  accessible  region.  Joint  efforts  to 
combine  simultaneous  observations  by  different  spacecraft 
and/or  ground-based  facilities  have  become  relatively  com¬ 
mon  recently.  The  problem  at  hand  is  how  to  utilise  the 
information  from  the  various  observations  to  reconstruct,  on 
a  global  scale,  the  instantaneous  distributions  of  all  relevant 
electrodynamical  parameters  in  a  self-consistent  marmet. 

Numerous  simulation  studies  of  the  interaction  between 
the  global  ionospheric  electric  field,  currents,  and  conductiv¬ 
ities  have  been  presented.  In  such  models,  two  of  the  quan¬ 
tities  are  specified,  and  the  third  is  calculated  from  these. 
The  chosen  input  quantities  are  typically  the  ionospheric 
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conductivity  combined  with  either  the  ionospheric  potential 
or  the  ionospheric  or  field-aligned  current  [e.g.,  Nisbet  et  eJ., 

1978;  Kamide  and  Matsushita,  1979a, b;  BJeuIer  et  al.,  1982; 
MarkJund  et  ai.,  1985].  The  conductivity  model  is  crucial  in 
models  of  this  kind  [cf.  Reiff,  1984]  and  can  therefore  not 
be  chosen  without  taking  into  account  the  physical  coupling 
between  the  conductivity  and  the  other  specified  input,  e.g., 
the  field-aligned  current.  This  has  usually  been  overlooked 
in  studies  of  this  kind  with  a  few  exceptions  [e.g.,  Mishin  et 
ai.,  1986] 

From  observations  it  is  clear  that  an  enhancement  in  the 
upward  field-aligned  current  is  usually  accompanied  by  an 
enhancement  in  the  ionisation  or  the  conductivity.  Although 
the  enhancements  of  these  quantities  are  typically  caused  by 
particles  belonging  to  different  energy  ranges,  such  a  rela¬ 
tionship  exists,  and  it  is  important  to  take  it  into  account 
in  the  modeling.  The  importance  of  the  coupling  between 
the  conductivity  and  the  upward  field-aligned  current  is 
treated  in  detail  in  a  special  report  [BJomberg  and  Mark- 
land,  1988s]. 

During  the  last  decade,  efforts  have  been  made  to  coor¬ 
dinate  observations  obtained  simultaneously  by  spacecraft 
and  ground-based  instrumentation  in  order  to  obtain  global 
snapshots  of  some  parameter,  often  the  electric  field  or  the 
ion  drift.  Schematic  convection  patterns  have  in  this  way 
been  constructed  for  the  nightside  high-latitude  ionosphere 
[HeeJis  et  al.,  1983]  and  for  the  dayside  ionosphere  [Mark- 
fund  et  ai.,  1986]. 

A  combination  of  this  latter  approach  (based  on  coordi¬ 
nation  of  simultaneous  observations)  with  that  using  global 
numerical  simulations  results  in  a  potentially  extremely  use¬ 
ful  way  to  better  understand  the  instantaneous  behavior  of 
the  magnetosphere-ionosphere  system.  This  is  so,  because 
not  only  one  single  parameter  (e.g.,  the  ion  convection  as 
exemplified  above)  but  all  electrodynamical  parameters  are 
in  a  self-consistent  way  included  in  such  an  approach.  In 
this  context  the  auroral  UV-imaging  instruments  found  on 
polar-orbiting  satellites  such  as  Dynamics  Explorer  I ,  Hilat, 
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and  Viking  play  a  crucial  role.  The  globed  view  of  the  aurora 
provides  a  natural  reference  frame  for  the  other  observations 
and  the  modeling,  and  it  reveals  important  information  on 
the  electrodynamical  state  of  the  different  parts  of  the  au¬ 
roral  ionosphere. 

Kamide  et  aj.  [1986]  used  images  of  auroral  emissions  ob¬ 
served  by  Dynamics  Explorer  1  to  calculate  high-latitude 
ionospheric  conductivities.  By  combining  these  with  simul¬ 
taneous  ground-based  magnetometer  data  they  were  able 
to  estimate  on  a  global  scale  the  distributions  of  ionospheric 
electric  fields  and  currents,  field-aligned  currents,  and  Joule 
heating  with  a  time  resolution  of  12  min. 

In  the  present  paper,  a  new  and  different  method  to  ob¬ 
tain  snapshots  of  the  auroral  electrodynamics  on  a  global 
scale  will  be  described,  using  auroral  UV  images  and  in  situ 
observations  of  fields  and  particles  (and/or  remote  observa¬ 
tions  of  the  same  or  related  parameters  by  ground-based 
instrumentation).  A  preliminary  and  brief  description  of 
this  new  technique  was  given  in  the  paper  by  Marklund  et 
aJ.  [1987]. 

In  the  present  paper  the  new  technique  and  numerical 
model  are  described  in  much  more  detail.  Furthermore,  a 
number  of  significant  improvements  and  additions  have  been 
made  in  the  analysis  of  the  particular  event  chosen  to  illus¬ 
trate  the  new  method.  The  most  important  ones  are  as 
follows; 

1.  We  included,  in  addition  to  the  Viking  data,  mag¬ 
netic  field  and  particle  data  from  a  noon-midnight  auroral 
oval  crossing  by  the  DMSP  F7  satellite.  As  a  result  of  this, 
all  information  used  to  calibrate  the  model  input  data  was 
obtained  within  15  min  for  the  dayside  oval  and  27  min  for 
the  entire  auroral  oval. 

2.  We  included  additional  data  from  the  other  experi¬ 
ments  on  Viking  (particle  and  high-frequency  wave  experi¬ 
ments)  which  made  possible  a  detailed  investigation  of  two 
auroral  acceleration  regions  (including  electric  held,  current, 
particle,  and  wave  signatures).  Furthermore,  a  more  de¬ 
tailed  estimate  of  the  downward  electron  energy  flux  at  the 
ionospheric  level  could  be  made. 

3.  We  calculated  the  E  x  B  drifts  from  the  Viking  elec¬ 
tric  field  data  which  allowed  a  two-dimensional  comparison 
with  the  model  convection  pattern. 

4.  We  calculated  the  ionospheric  currents.  Joule  power 
dissipation,  and  convection  in  the  equatorial  plane  of  the 
magnetosphere. 

Section  2  outlines  the  main  principles  of  how  to  construct 
the  model  input  data  from  the  observations.  Section  3  de¬ 
scribes  the  numerical  model  (in  particular  the  choice  of  con¬ 
ductivity  model).  The  auroral  situation,  the  geophysical 
conditions,  and  the  geometry  of  the  oval  are  described  in 
section  4.  Section  5  presents  the  observations,  model  input 
distributions,  and  results  of  the  model  calculations.  Discus¬ 
sion  and  summary  of  the  results  are  presented  in  sections  6 
and  7,  respectively. 

2.  Methodology 

The  basic  elements  needed  to  infer  and  calibrate  the  model 
input  data  and  check  the  output  data,  following  the  proce¬ 
dure  used  heic,  are  (1)  one  or  several  UV  images  showing 
clearly  the  auroral  emissions  along  the  entire  auroral  oval 
as  well  as  in  the  polar  cap,  (2)  in  situ  observations  of  elec¬ 
tric  and  magnetic  fields  and  auroral  particles  from  polar- 
orbiting  satellites  (at  least  one),  preferably  supplemented 


by  remote  observations  of  these  or  related  parameters  by 
ground-based  facilities  such  as  incoherent  or  coherent  scat¬ 
ter  radars  and  magnetometer  networks,  (3)  experience  of 
the  statistical  patterns  of  field-aligned  currents  and  iono¬ 
spheric  conductivities  and  the  observed  coupling  between 
these  for  different  geophysical  conditions,  as  reported  in  the 
literature,  and  (4)  a  sophisticated  numerical  model  to  solve 
the  elliptical  equation  that  relates  the  ionospheric  potential 
(model  output)  to  the  ionospheric  conductivity  and  field- 
aligned  current  distributions  (model  input). 

The  task  is  to  reconstruct  on  a  global  scale  the  distribu¬ 
tions  of  all  relevant  electrodynamic  parameters  associated 
with  the  particular  auroral  situation  given  by  the  UV  image. 
To  accomplish  this,  global  input  distributions  of  conductiv¬ 
ities  and  magnetic  field-aligned  currents  have  to  be  inferred 
from  and  made  consistent  with  the  direct  observations  and 
fed  into  the  numerical  model  to  obtain  the  potential,  which 
is  finally  checked  against  the  observations. 

The  Auroral  UV  Image 

The  auroral  distribution  viewed  by  the  UV  imager  is  used 
to  obtain  a  global  reference  baine  for  the  other  observations 
and  in  particular  to  locate  the  regions  of  active  intense  au¬ 
roras.  These  are  clearly  associated  with  enhanced  upward 
field-aligned  currents  and  ionospheric  conductivities.  The 
coupling  between  these  quantities  is  represented  here  by  the 
linear  relationship  expressed  by  equation  (10)  below  (in  sec¬ 
tion  3). 

The  UV  images  also  reveal  much  fine  structure  both  in 
the  intense  auroral  structures  and  in  the  weaker  and  diffuse 
“background”  auroral  oval.  The  fine  structure  in  the  auro¬ 
ral  emissions  has  not  been  taken  into  account  in  the  global 
modeling  for  several  obvious  reasons.  (1)  There  are  several 
difficulties  involved  in  trying  to  obtain  reliable  absolute  in¬ 
tensities  from  the  UV  image.  For  example,  corrections  have 
to  be  made  for  effects  of  the  oblique  viewing  angle  of  the 
imager  and  for  different  contrasts  in  different  parts  of  the 
auroral  oval.  (2)  The  errors  arising  from  trying  to  take  the 
fine  structure  of  the  diffuse  background  auroral  oval  into  ac¬ 
count  in  the  modeling  would  most  likely  be  comparable  to 
(or  even  larger  than)  the  real  intensity  variations.  (3)  Even 
if  it  would  be  possible,  in  principle,  to  infer  some  fine  struc¬ 
ture  in  the  ionospheric  conductivity  from  the  UV  image, 
this  would  not  improve  the  modeling  unless  the  field-aligned 
current  distribution  could  be  inferred  with  a  corresponding 
spatial  resolution  (which  is  clearly  not  possible). 

We  have  therefore  chosen  to  represent  the  diffuse  “back¬ 
ground”  auroral  oval  by  statistical  patterns  of  conductivi¬ 
ties  and  field-aligned  currents  .  The  MLT  dependence  of 
this  field-aligned  current  distribution  is  taken  from  Jijima 
and  Potemra  [1976a,  J978],  while  its  location  is  chosen  in 
consistency  with  the  UV  image. 

To  summarise,  the  UV  image  is  used  to  obtain  informa¬ 
tion  on  the  sise,  shape,  and  width  of  the  oval  as  a  general 
frame  of  reference  for  the  other  observations  and  for  the 
modeling.  Moreover,  the  UV  image  is  used  to  identify  the 
location  and  extent  in  latitude  and  local  time  of  the  more 
intense  auroral  emissions  along  the  oval.  Statistical  pat¬ 
terns  of  field-cdigned  current  represent  the  current  system 
associated  with  the  diffuse  aurora,  and  the  corresponding 
conductivity  distribution  is  given  by  the  Gaussian-shaped, 
“oval  centered,”  term  medified  by  enhancements  in  the  re¬ 
gions  of  statistical  upward  field-aligned  currents.  In  addi- 
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tion  to  this,  the  current  and  conductivity  in  regions  where 
active  aurora  is  found  from  the  image  are  enhanced.  Note 
that  for  both  these  contributions,  the  current  and  conduc¬ 
tivity  are  treated  in  a  consistent  manner  in  the  sense  that 
part  of  the  conductivity  distribution  is  coupled  to  the  field- 
aligned  current  (cf.  section  3).  Moreover,  the  field-aligned 
current  and  the  conductivity  distribution  are  chosen  so  as  to 
be  everywhere  continuous.  The  center  latitude  between  the 
region  1  and  region  2  field-aligned  currents  and  that  of  the 
peak  background  conductivity  coincide  in  the  model  with 
the  center  latitude  of  the  auroral  oval  as  observed  by  the 
UV  imager. 

Calibration  of  the  Model  Input 

The  rough  and  qualitative  picture  of  the  model  input  dis¬ 
tributions  resulting  from  the  procedure  above  will  now  be 
refined  and  calibrated  to  fit  with  the  actual  measurements 
(in  situ  and/or  remote).  For  the  sake  of  convenience,  the 
calibration  procedure  will  be  described  here  using  satellite 
observations,  although  it  applies  equally  well  to  ground- 
based  observations.  Magnetic  field  and  particle  data  are 
used  to  calculate  the  distributions  of  field-aligned  currents 
and  particle  energy  flux  along  the  satellite  orbit.  These  es¬ 
timates  have  to  be  multiplied  with  the  appropriate  scaling 
factors  between  the  satellite  altitude  and  the  ionosphere.  In 
this  study  we  have  used  the  scaling  factors  obtained  from  a 
pure  dipole  magnetic  field. 

The  particle  energy  flux  is  used  to  calculate  the  contribu¬ 
tion  from  precipitating  particles  to  the  ionospheric  conduc¬ 
tivity  as  expressed  by  the  formula  Sp  =  from  Harel 

et  al.  [1981]  where  o  is  a  proportionality  constant  and  F 
the  particle  energy  flux.  Another,  perhaps  more  accurate 
way  to  calculate  Dp  is  to  use  the  dissipation  algorithm  for 
electrons,  given  by  Rees  [1963]  assuming  a  balance  between 
ionization  and  recombination  losses.  Since  Sp  does  not  rep¬ 
resent  the  total  conductivity  but  only  the  contribution  due 
to  the  precipitating  particles,  it  should  be  used  to  calibrate 
the  corresponding  part  of  the  model  conductivity  (i.e.,  ex¬ 
cluding  the  contribution  to  the  conductivity  from  solar  EUV 
radiation;  see  section  3). 

For  orbits  or  segments  of  orbits  for  which  there  is  a  good 
large-scale  correlation  between  orthogonal  components  of 
the  electric  field  (E)  and  the  residual  magnetic  field,  b,  the 
total  ionospheric  conductivity,  Sp,  can  be  inferred  [Sug- 
iura,  1984]  according  to  the  following  expression:  Sp  = 
Ab/^oAE.  This  constitutes  another  possible  wav  to  cali¬ 
brate  the  model  input  conductivity. 

Calibration  of  the  model  input  data  is  now  possible  us¬ 
ing  these  estimates  of  the  field-aligned  current  and  height- 
integrated  conductivity  distributions.  Note  that  all  model 
coefficients  related  to  the  background  conductivity  (cf.  sec¬ 
tion  3)  are  position-independent.  Thus,  the  calibration  of 
the  background  conductivity  is,  with  necessity,  made  glob¬ 
ally. 

For  regions  covered  by  the  satellite  orbit  the  model  cur¬ 
rents  and  conductivities  are  matched  to  the  observed  widths 
of  the  large-scale  current  sheets  and  particle  precipitation 
regions  and  to  the  typical  magnitudes  of  these  parameters. 
No  attempt  has  been  made  to  represent  the  observed  fine 
structure  in  the  model  input  data. 

For  active  regions  not  covered  by  the  satellite  orbit(s), 
both  the  region  1  and  region  2  large-scale  model  currents 
are  enhanced  so  that  their  amplitudes  are  consistent  with 


either  (1)  current  magnitudes  measured  in  some  nearby  ac¬ 
tive  auroral  structure  with  a  similar  UV  emission  inten¬ 
sity  or  (2)  typical  current  magnitudes,  based  on  the  experi¬ 
ence  from  earlier  satellite  measurements  over  similar  auroral 
structures. 

The  region  2  to  region  1  current  ratio  will  in  this  way 
assume  a  lower  value  in  the  regions  of  low  activity  (namely 
the  ratio  specified  by  the  quiet  time  statistical  patterns  of 
lijima  and  Poteirua  [1976a,  1978]  )  than  in  the  regions  of 
active  aurora  where  this  ratio  approaches  unity. 

Since  the  field-aligned  current  and  the  conductivity  are 
coupled  in  the  model,  as  represented  by  one  of  the  terms 
in  the  expression  for  the  conductivity  (cf.  section  3),  the 
calibration  of  the  conductivity  becomes  an  iterative  process 
during  which  a  better  estimate  of  the  proportionality  factor 
in  equation  (10)  will  evolve.  A  complete  description  of  the 
conductivity  model  is  given  in  section  3. 

Model  Output  and  Check  Procedure 

Once  the  model  input  data  are  specified,  they  are  fed  into 
the  numerical  model  (described  in  section  3)  which  calcu¬ 
lates  the  global  ionospheric  potential  distribution.  Before 
presenting  the  potential  pattern  as  it  comes  out  from  the 
calculation  the  result  is  to  be  checked  agsunst  measurements 
in  the  same  way  as  the  model  input  data. 

Electric  field  (or  ion  drift)  observations  are  used  to  calcu¬ 
late  both  the  potential  variation  and  the  E  x  B  drift  pattern 
along  the  projected  satellite  orbits  (or  radar  tracks),  which 
allows  point-by-point  comparisons  with  the  model  result. 

Since  the  strengths  of  the  field-aligned  currents  (and  the 
corresponding  conductivity)  in  active  auroral  regions  not 
covered  by  the  satellite  are  assumed  rather  than  measured, 
they  could  be  varied  in  such  a  way  that  a  better  fit  is  ob¬ 
tained  between  the  measured  and  computed  potential  vari¬ 
ation.  Note,  however,  that  such  a  variation  in  the  current 
strength  of,  say,  a  nightside  amoral  structure  will  have  only 
minor  influence  on  the  dayside  potential  pattern  provided 
that  not  too  extreme  current  values  have  been  used.  More 
locally,  the  effects  are,  however,  significant  and  could  be 
used  for  improvement  of  the  model  input  data. 

The  resulting  global  pattern  can  be  used  to  predict  the 
electric  field  or  convection  velocity  along  the  ionospheric 
foot  point  of  any  satellite  trajectory  or  radar  track  which  en¬ 
ables  direct  point-by-point  comparisons  with  observations. 
This  should  be  extremely  useful  in,  for  example,  coordi¬ 
nated  multiradar-multisatellite  studies  of  the  global  con¬ 
vection.  In  such  a  case,  with  a  relatively  dense  network 
of  observations,  the  iterative  procedure,  to  improve  succes¬ 
sively  the  fits  between  the  measured  and  computed  potential 
variations  along  the  different  satellite  trajectories  (or  radar 
tracks),  as  described  above,  should  be  most  valuable. 

Combining  the  electric  field  given  by  the  final  potential 
pattern  with  the  model  conductivity  allows  both  the  iono¬ 
spheric  currents  and  the  Joule  power  dissipation  to  be  calcu¬ 
lated  according  to  equations  (I)  and  (11)  below.  The  effect 
of  these  calculated  ionospheric  currents  at  the  ground  level 
can  be  estimated,  to  be  used  for  direct  comparisons  with 
data  from  ground-based  magnetometers.  This  represents 
another  possible  way  to  check  the  model  data. 

An  important  feature  of  the  new  technique  presented  here 
is  that  all  relevant  electrodynamical  parameters  involved  in 
the  modeling,  namely  the  electric  field,  the  field-aligned  cur¬ 
rents,  and  the  conductivities,  are  treated  self-consistently 
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and  checked  ot  calibiated  point  by  point  against  direct  mea¬ 
surements  in  an  interactive  way  not  done  before. 


(To) 


3.  Model  Equations 

The  calculation  of  the  potential  distribution  has  been  per¬ 
formed  using  spherical  (surface)  coordinates,  where  the  9 
coordinate  corresponds  to  magnetic  colatitude  and  the  ^ 
coordinate  is  related  to  magnetic  local  time  (MLT). 

The  basic  equation  relating  the  height-integrated  iono¬ 
spheric  (horizontal)  current  (Jo)  to  the  ionospheric  electric 
field  (Eo)  can  be  written 


Jo  —  E  •  (Eo  +  VnXB)  (1) 


where  Eo  is  the  electric  field  in  the  Earth-fixed  frame  and 
v„  is  the  velocity  of  the  neutral  wind,  which  in  the  calcu¬ 
lation  below  has  been  assumed  to  be  zero.  Note  that  the 
electromotive  force  Vn  x  B  does  not  significantly  influence 
the  field  -aligned  current  [cf.  Lyons  and  WaJteracheid,  1986]. 

E  is  the  height-integrated  horizontal  conductivity  tensor 


E  = 


Ess 

Ess 


(2) 


Its  components  (except  at  equatorial  latitudes)  are  given  by 


Sh  =  (76) 

The  root-sum-square  relation  is  employed,  since  the  vari¬ 
ous  ionization  processes  giving  rise  to  the  conductivity  take 
place  roughly  within  the  same  altitude  range,  and  hence, 
as  discussed  by,  for  example,  Wallis  and  Budzinski  [1981], 
addition  of  the  ionisation  rates  (which  if  balzuiced  by  recom¬ 
bination  are  proportional  to  n^)  gives  a  far  more  accurate 
estimate  of  the  total  conductivity  (proportional  to  n«)  than 
a  direct  summation  of  the  different  conductivity  contribu¬ 
tions. 

So  is  a  background  term  having  a  constant  low  value,  rep¬ 
resenting  the  contribution  from  cosmic  radiation  and  galac¬ 
tic  EUV.  The  exact  magnitude  of  this  term  is  not  critical  to 
the  modeling  and  has  in  this  study  been  given  the  value  0.5 
S. 

Ef  represents  an  average  auroral  zone  background  con¬ 
ductivity  given  by 


E,  =  E,o  exp 


( 


(9  -  9oy 


) 


(8) 


Ess  =  Sp/sin^  /  (Ju) 

Ess  “  —Ess  =  Eff  / sin  1  ( J^) 

Ess  =  Sp  (3c) 

where  I  is  the  inclination  of  the  geomagnetic  field  (B)  (ap¬ 
proximated  by  a  dipole)  and  Ep  and  Eh  are  the  height- 
integrated  Pedersen  and  Hall  conductivities,  respectively. 

The  height-integrated  horizontal  current  and  the  field- 
aligned  current  (positive  downward)  are  related  through 


V  •  Jo  =  ill  •  sin/  (4) 

V  (EV*)= -ill -sin/  (5) 

This  equation  forms  the  basis  for  solving  (numerically)  the 
electrostatic  potential  once  the  field-aligned  current  and  the 
conductivity  are  specified. 

The  boundary  condition  used  in  the  calculation  is  #(9  = 
90”,  1^)  =  0;  i.e.,  the  magnetic  equator  is  an  equipotential 
surface. 

In  the  present  model  the  ordinary  region  1  and  region  2 
field-aligned  currents  [cf.  lijima  and  Potemra,  19768,  1978], 
as  well  as  the  cusp  currents  [Sjima  and  Potemra,  1976  b]  can 
be  represented.  The  basic  parameters  defining  the  geome¬ 
try  of  the  current  systems  are  ^oo>  colatitude  of  the  dayside 
region  1  /  region  2  interface;  ^9o,  difference  between  night- 
side  and  dayside  interface  latitudes;  and  A^i,  Affc, 

the  latitudinal  widths  (in  degrees)  of  region  1,  region  2,  and 
cusp  currents,  respectively. 

The  center  of  the  auroral  oval  (assumed  to  coincide  with 
the  region  1  /  region  2  interface)  is  described  in  this  study 
by 


9o  —  9oo  ■+■  * 


cos 


15”  •  MLT 
2 


(6) 


The  height-integrated  conductivities  are  represented  by 


centered  around  the  colatitude  9o  (coinciding  with  the  center 
of  the  auroral  oval)  and  having  a  characteristic  half  width 
of  9f0  degrees. 

Ei/v  represents  the  part  of  the  conductivities  which  is 
produced  by  solar  EUV  radiation.  It  is  given  here  by 

Edv  =  Ofy/cos  X  (9) 

where  x  fbe  solar  zenith  angle  and  a  a  proportionality 
factor.  The  value  of  a  is  taken  from  Vickrey  et  al.  [1981]. 

Ej|l  represents  the  conductivity  enhancement  produced  by 
precipitating  electrons  associated  with  upward  field-aligned 
currents.  In  the  present  model  a  relationship  between  the 
field-aligned  current  and  the  conductivity  has  been  used, 
such  that  the  conductivity  peaks  in  regions  of  upward  cur¬ 
rent.  The  form  of  this  relation  is 

_  f  fc(MLT)|j|||  j||  upward 
'N  ( 0  y||  downward  '  ' 

The  factor  of  proportionality  k(MLT)  is  here  chosen  to  be 
dependent  of  local  time  in  order  to  account  for  the  differ¬ 
ences  in  the  hardness  of  the  particle  spectrum.  Thus,  it  has 
a  lower  value  on  the  dayside,  where  the  precipitation  gen¬ 
erally  is  softer  than  on  the  nightside.  In  the  morning  and 
evening  sectors  there  is  a  smooth  transition  between  the  two 
extreme  values.  The  transition  region  stretches  over  2  hours 
of  magnetic  local  time.  The  contribution  from  protons  to 
the  total  precipitated  particle  energy  is  estimated  to  be  at 
most  10  to  20  per  cent  [Huftqvist,  1973]  and  is  assumed 
here  to  be  accounted  for  by  E^.  An  important  parameter 
in  the  modeling  is  the  Hall  to  Pedersen  conductivity  ratio 
Eh/Ef>  which,  in  principle,  should  be  slightly  different  for 
the  various  terms  in  the  conductivity  model  (cf.  equation 
(7))  due  to  different  ionization  mechanisms.  The  ratio  be¬ 
tween  the  total  Eh  and  Ep  used  in  this  study  is  1.7  within 
the  auroral  oval  and  close  to  1.0  in  the  nightside  polar  cap 
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TABLE  1.  Model  Input  Parameter  Values 


Geometry  of  the  Model  Currents 


#00 

=  H” 

A#o 

=  12° 

=  Aflj  =  3° 

A/lc 

=  2° 

Conductivity  Model 

Lop 

=  0.5  S 

Lgop 

=  5  S 

#0 

=  14°  at  1200  MLT 

9(0 

=  3®  (cf.  equation  (8)) 

ap 

=  5  S  (cf.  equation  (9)) 

f  2  S/(/iA/m*),  on  dayside 

#ep 

\  13  S/(/iA/m^),  on  nightside 

Lh/Lp 

f  1,  for  Eo 

“  ^  1.7,  for  SuVi  *nd 

Corotation  Potential 

P 

=  8.1  X  10^’  Am*  (cf.  equation  (12)) 

Ul 

=  2ir/(24  X  3600) 

=  11°  (cf.  equation  ^13)) 

#NP 

UT 

=  0800 

For  explanation,  see  section  3. 


and  subauroral  regions  using  the  parameter  values  listed  in 
Table  1. 

Once  the  ionospheric  potential  is  calculated  from  equation 
(5),  the  associated  electric  field  and  the  height-integrated 
conductivity  can  be  used  to  calculate  the  ionospheric  cur¬ 
rent,  Jo ,  given  by  ( 1 )  and  the  Joule  power  dissipation, 
given  below: 

=  B  e  e  =  +  (11) 

The  ionospheric  potential,  4,  can  be  projected  onto  the  mag- 
netospheric  equatorial  plane  under  the  assumption  of  perfect 
mapping  (no  magnetic  field-aligned  electric  fields).  In  the 
present  study  a  simple  centered  dipole  field  has  been  used. 

The  electrostatic  potential  or  the  electric  field  can  be  given 
in  either  the  corotating  frame  of  reference  (assumed  to  be 
identical  to  the  frame  where  Vn  =  0)  or  in  a  coordinate  sys¬ 
tem  defined  by  the  magnetic  axis  and  the  direction  to  the 
Sun  (Sun-fixed  frame).  The  geographic  north  pole  can  be 
regarded  as  rotating  around  the  magnetic  pole  at  a  mag¬ 
netic  colatitude  of  11°  with  a  24-hour  period.  The  corota¬ 
tion  potential  (or  corotation  electric  field)  needed  for  trans¬ 
formation  between  the  two  frames  will  therefore  be  UT- 
dependent. 

To  account  for  this  in  the  model  an  expression  has  been 
derived  for  the  UT-dependent  corotation  electric  field,  which 
in  turn  has  been  used  to  obtain  a  “pseudo-potential,”  ♦c, 
which  approximates  the  effect  of  the  UT-dependence  of  the 
corotation  field.  This  “pseudo-potential”  reads: 

♦c(®,  ^)  =  sin0Np  cos(^NP  -  ^)  (^  -(-  —z — ) 

—  cos^NP  sin*  S  (12) 

where  p  is  the  magnetic  dipole  moment  (=  8.1  x  10**  Am*); 
it>  is  the  angular  velocity  of  the  Earth’s  rotation;  r  is  the 


geocentric  distance;  and  and  are  magnetic  colatitude 
and  local  time,  respectively,  with  subscript  NP  referring  to 
the  geographic  north  pole.  The  latter  are  approximately 
given  by 

^NP  «  11°  (13o) 

I^NP  =  15  •  UT-I- 111°  (136) 

An  exhaustive  treatment  of  the  corotation  field  and  its  UT- 
dependence  will  be  the  subject  of  a  future  study.  For  a 
further  description  of  the  numerical  method,  see  Blomberg 
and  Marklund  [1988b].  The  values  of  the  model  parameters 
used  for  this  study  are  summarised  in  Table  1  . 

4.  Geophysical  Conditions 

Figure  1  shows  the  northern  hemisphere  auroral  oval  as 
viewed  by  the  UV  imager  experiment  on  Viking  at  0T24 
UT  during  orbit  360,  April  28,  1986.  The  auroral  oved 
was  crossed  by  Viking  in  roughly  a  dusk-dawn  orbit  and 
by  DMSP  FT  in  roughly  a  noon-midnight  orbit,  as  shown 
by  the  inserted  ionospheric  foot  points  of  the  two  satellites. 

Three  localised  bright  nuroral  structures  can  be  recog¬ 
nised  in  the  evening  and  midnight  sector,  the  largest  one 
centered  at  about  65°  invariant  latitude,  0100  magnetic  lo¬ 
cal  time  (MLT)  with  an  east- west  extent  of  roughly  20°. 
The  western  edge  of  the  bright  spot  around  2130  MLT  was 
crossed  by  the  DMSP  FT  satellite,  as  will  be  shown  below. 
Relatively  bright  dayside  auroras  are  seen  in  the  ptenoon 
sector,  whereas  the  remaining  part  of  the  auroral  oval  has  a 
relatively  smooth  and  regular  luminosity  distribution.  Note 
that  the  auroral  oval  in  the  early  morning  sector  nearly  co¬ 
incides  with  the  limb  of  the  Earth.  This  implies  that  the 
weak  emission  intensities  observed  in  this  region  may  not  be 
due  entirely  to  auroral  precipitation. 

The  afternoon  sector  of  the  auroral  oval  was  traversed  by 
Viking  around  0T24  UT,  i.e.,  the  time  the  UV  image  was 
taken.  The  prenoon  sector  of  the  auroral  oval  was  traversed 
by  the  DMSP  FT  satellite  around  0T39  UT;  i.e.,  all  informa¬ 
tion  used  to  reconstruct  the  dayside  model  input  data  was 
obtained  roughly  within  a  15-min  interval. 

Of  the  nightside  auroral  activity,  the  bright  spot  centered 
around  2130  MLT  was  traversed  by  the  DMSP  FT  satellite 
around  0T50  UT,  which  allowed  a  calibration  of  the  model 
input  data  for  this  part  of  the  nightside  oval  within  less  than 
half  an  hour  after  the  UV  image  was  taken. 

The  AE  index  for  this  time  period  [National  Space  Science 
Data  Center,  1986]  shows  that  after  a  period  of  growth  in  the 
electrojet  intensity  previous  to  the  UV  image,  it  remained  at 
a  high  but  relatively  stable  value  around  500  nT  during  the 
calibration  time  period  0T24  -  0T51  UT  and  then  declined 
to  about  400  nT  as  Viking  traversed  the  prenoon  auroral 
oval  around  0900  UT. 

We  have  chosen  this  particular  event  to  illustrate  in  some 
detail  how  to  use  the  new  technique  to  reconstruct  the  elec¬ 
trodynamics  associated  with  a  given  auroral  situation  viewed 
by  a  UV  imager.  The  assumption  made  to  enable  such 
a  reconstruction  on  a  global  scale  is  that  the  auroral  and 
geophysical  conditions  remained  relatively  stable  during  the 
time  period  used  for  the  calibration  of  the  model  input  data, 
i.e.,  between  OT24  UT  and  0T50  UT  (for  the  dayside  auroral 
oval  between  0T24  UT  and  0T39  UT).  Since  in  situ  observa¬ 
tions  of  the  nightside  auroral  activity  are  available  only  from 
one  oval  crossing  (DMSP  FT,  2130  MLT,  OT39  UT)  as  com- 
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Fig  1.  The  northern  hemisphere  Auroral  oval  as  viewed  by  the 
UV  imager  experiment  on  Viking  at  0724  UT,  orbit  360,  April 
28,  1986 


pared  to  three  on  the  dayside,  the  possibility  of  calibrating 
the  input  data  and  checking  the  stability  is  smaller  there. 
The  DMSP  F7  observations  of  strongly  enhanced  particle 
precipitation  and  field -aligned  currents  (which  will  be  pre¬ 
sented  below)  at  the  same  location  at  which  the  intense  and 
relatively  narrow  auroral  structure  around  2130  MLT  was 
observed  by  the  UV  imager  27  min  earlier  speak  in  favor 
of  a  less  violent  nightside  auroral  activity  during  thb  time 
period. 

No  IMP  8  interplanetary  magnetic  field  (IMF)  data  exist 
for  this  time  period,  but  data  from  the  Greenland  magne¬ 
tometer  chain  indicate  that  the  IMF  B,  component  was  pos¬ 
itive  during  this  event  (E.  Friis-Christensen,  private  coiiunu- 
nication,  1987). 

5.  Results  and  Comparisons  With  Viking  and 
DMSP  F7  Satellite  Data 

Figure  2  shows  the  input  Sp  and  y||  distributions,  which 
are  qualitatively  consistent  with  the  auroral  image  seen  in 
Figure  1  and,  as  will  be  shown  below,  in  good  quantitative 
agreement  with  the  Viking  and  DMSP  F7  data. 

The  region  1  and  region  2  current  patterns  of  lijima  and 
Potemra  [1976a,  1978],  representative  for  the  diffuse  quiet 
time  auroral  oval,  are  seen  to  be  modified  by  (1)  current 
enhancements  in  the  local  time  sectors  associated  with  the 
bright  auroral  structures  identified  from  Figure  1  and  (2) 
cusp  current  sheets,  inferred  from  the  Viking  magnetic  field 
observations  (see  below)  with  location  and  extent  in  accor¬ 
dance  with  the  observations  by  lijima  and  Potemra  [1976b]. 
Since  both  the  upward  and  downward  currents  associated 


with  a  bright  auroral  structure  are  assumed  to  be  enhanced, 
the  ;(|  contours  will  look  relatively  symmetric  about  the 
noon-midnight  plane  in  contrast  to  the  Dp  distribution, 
which  is  seen  to  have  a  relatively  higher  magnitude  in  the 
regions  of  upward  field-aligned  currents  (as  a  result  of  the 
current-conductivity  coupling  term)  with  intense  maxima  in 
the  bright  auroral  forms.  The  ionisation  contribution  from 
solar  EUV  gives  rise  to  the  day-night  conductivity  gradient. 

The  two  projected  orbits  are  seen  to  cross  the  prenoon 
auroral  oval  (region  1  and  region  2)  relatively  close  and  to 
cross  each  other  almost  at  the  center  of  the  prenoon  cusp 
current  sheet.  The  1-hour  time  separation  between  the  two 
satellites  at  this  location  allows  us  to  check  the  time  stability 
of  this  prenoon  triple  current  sheet  (see  below).  Note  that 
the  DMSP  F7  satellite  passed  through  the  intense  bright 
auroral  form  around  2130  MLT. 


12 
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Fig.  2.  Input  Pedersen  conductivity  (Ey>,  top)  and  input  field- 
aligned  current  (j||,  bottom).  Included  are  also  the  magnetic  foot 
points  of  the  Viking  (dusk-dawn)  and  DMSP  F7  (noon-midnight) 
satellites. 
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Pig.  3.  Quantitative  point-by-point  comparison  between  model  input  data  and  Viking  data  along  the  ionospheric 
projection  of  the  Viking  orbit. 

Figure  3  shows,  at  the  top,  the  model  and  calculated  portional  to  the  square  root  of  the  downward  electron  en- 
height-integrated  Pedersen  conductivity  profiles  and,  at  the  ergy  flux  [Harel  et  of,  1981]  as  described  in  section  2.  The 
bottom,  the  model  and  calculated  field-aligned  current  pro-  fleid-aligned  current,  y||,  has  been  estimated  &om  the  resid- 
files  along  the  ionospheric  projection  of  the  Viking  orbit,  ual  of  the  spin  axis  component  of  the  magnetic  field,  6j, 
The  contribution  to  Ep  from  the  measured  particles  has  shown  in  Figure  4  (bottom)  and  obtained  by  subtraction  of 
been  estimated  using  the  assumption  that  it  is  roughly  pro-  the  slow  time  variation  of  the  background  magnetic  field. 
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Fig.  4.  Electric  field  (top)  and  magnetic  field  data  (bottom)  obtained  for  the  dusk-dawn  Viking  orbit  360.  For 
a  definition  of  the  coordinate  system  used,  see  Block  et  ml.  (1987],  Ej  corresponds  roughly  to  the  dawn-dusk 
electric  field,  and  Bs  to  the  sunward  magnetic  field. 
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Fig.  5.  Quantitative  point -by-point  comparison  between  model  input  data  and  OMSP  P7  data  along  the  iono¬ 
spheric  projection  of  the  DMSP  F7  orbit. 


(This  component  contained  the  major  part  of  the  smatler- 
scale  variations  seen  in  the  Viking  magnetic  field  data.)  For 
a  definition  of  the  coordinate  system  used,  see  Block  et  at. 
[1987],  The  current  values  have  been  scaled  using  a  dipole 
magnetic  field  approximation  so  that  they  refer  to  the  iono¬ 
spheric  level.  Note  the  clear  large-scale  correlation  between 
the  orthogonal  electric  field  and  magnetic  field  components 
in  Figure  4  .  This  allows  an  independent  estimate  of  the 
variation  (denoted  by  circles  in  Figure  3  )  as  demonstrated 
by  Sugiura  [1984). 

As  can  be  seen,  the  Up  and  jg  profiles  calculated  from  the 
Viking  data  are  very  well  represented  by  the  model  data  in 
terms  of  the  relative  magnitudes,  widths,  and  locations  of 
the  major  structures.  Note  that  the  calculated  conductivity 
always  has  local  maxima  in  the  regions  of  upward  field- 
aligned  currents.  This  field-aligned  current-conductivity 
coupling  is  a  characteristic  feature  in  auroral  electrodynam¬ 
ics  and  has  been  taken  into  account  here  by  using  the  linear 
relationship  between  jn  and  S,  described  by  equation  (10). 

Figure  5  shows  in  a  fashion  similar  to  that  in  Figure  3 
comparisons  between  the  model  field-aligned  current  and 
conductivity  profiles  and  the  corresponding  estimates  based 
on  the  DMSP  F7  data.  Here,  too,  the  agreement  is  good 
concerning  the  large-scale  variations.  Consider  for  exstm- 
ple  the  triple  current  sheet  encountered  around  0740  UT. 
The  location  and  the  integrated  values  of  the  upflowing  and 
downflowing  currents  are  very  similar,  although  the  width  of 
the  model  current  is  slightly  greater.  As  mentioned  above, 
this  triple  current  sheet  was  also  traversed  by  the  Viking 
satellite  about  1  hour  later.  The  magnitudes  and  signatures 
of  the  currents  remain  essentially  the  same  during  this  1- 
hour  time  period,  with  a  value  of  0.8  -  1  pA/m*  for  the  re¬ 
gion  2  upward  current,  3  for  the  region  1  downward 

current,  and  about  2  fik/vc?  for  the  upward  cusp  current. 
Thus,  u  far  as  one  could  judge  from  these  observations,  the 


dayside  conditions  seem  to  have  been  relatively  stable  even 
for  the  1-hour  time  period  preceding  the  DMSP  F7  crossing. 

For  the  nightside  auroral  crossing,  the  model  conductivity 
and  field-aligned  current  profiles,  representing  the  intense 
auroral  structure  around  2130  MLT  seen  in  the  UV  image, 
agree  relatively  well  with  the  corresponding  estimates  based 
on  the  DMSP  F7  data.  The  narrow  structure  at  the  pole- 
ward  edge  of  the  oval  (0750  UT)  is,  however,  not  represented 
in  the  model.  This  may  be  just  a  local  arc  structure,  or  it 
may  represent  the  poleward  part  of  a  multiple  current  sheet 
associated  with  the  Harang  discontinuity.  Adopting  the  for¬ 
mer  interpretation,  there  is  no  reason  to  include  small-scale 
structures  in  the  model  only  along  the  satellite  orbits  where 
these  can  be  detected  and  not  ebewhere.  The  net  current 
from  these  opposite  currents  is  close  to  lero  and  will  there¬ 
fore  not  affect  the  results  except  very  locally. 

To  summarue,  the  model  input  data  are  believed  to  rather 
accurately  represent  the  real  field-aligned  current  and  con¬ 
ductivity  dutributions,  as  indicated  by  the  good  agreement 
with  the  measurements  along  the  two  different  satellite  or¬ 
bits. 

The  equipotential  distribution,  calculated  using  the  nu¬ 
merical  model  described  in  section  3  (equation  (5))  and  the 
input  distributions  presented  above,  is  presented  in  Figure  6 
in  an  inertial  frame  of  reference.  Superposed  on  this  we 
have  for  comparbon  plotted  the  E  y  B  drift  vectors  calcu¬ 
lated  from  the  Viking  electric  field  (Figure  4  )  assuming  that 
E  ■  B  =  0.  Note  that  the  drift  vectors  are  almost  parallel  to 
the  convection  streamlines  everywhere  except  for  two  nar¬ 
row  regions,  as  will  be  discussed  below. 

The  calculated  convection  pattern  b  of  the  normal  two¬ 
cell  type  with  a  weak  dusk  cell  displaced  toward  noon  and  a 
more  intense,  longitudinally  extended  and  crescent-shaped 
dawn  cell  characteristic  for  the  prevailing  positive  IMF  By . 
The  pattern  is  quite  similar  to  the  BC  model  of  Heppner 
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Fig  6  Output  pot^nlial  distribution  in  the  inertial  frAm^  Con¬ 
tour  separation  is  5  kV  Along  the  Viking  foot  point  (dusk-dawn) 
are  shown  the  E  x  B  drift  vectors  calculated  from  Viking  electric 
field  measurements 

and  Afaynard  [1987],  representing  the  most  typical  potential 
pattern  encountered  in  the  northern  hemisphere  for  positive 
IMF  By  conditions  and  Kp  in  the  range  <  Kp  <  4' 
(samf  as  in  th?  present  study).  The  main  difference  is  that 
the  equipotential  contours  of  the  statistical  BC  model  ate 
parallel  to  the  0900-2100  MLT  meridian  within  almost  the 
entire  polar  cap  whereas  the  potential  contours  of  the  “in¬ 
stantaneous"  pattern  shown  in  Figure  6  change  from  being 
almost  Sun-aligned  on  the  dayside  polar  cap  to  aligned  in 
the  dawn-dusk  direction  on  the  nightside  polar  cap.  An¬ 


other  difference  is  that  the  morning  cell  dominates  more 
clearly  over  the  evening  cell  for  the  “instantaneous”  pattern. 

Note  the  westward  intrusion  or  clockwise  rotation  of  the 
convection  pattern  close  to  midnight.  Such  a  clockwise  rota¬ 
tion  is  caused  by  the  combined  effect  of  intense  field-aligned 
currents  and  polarisation  electric  fields  associated  with  the 
bright  auroral  intensifications  there  [Marklund  et  al.,  1985] 
as  well  as  the  day-night  conductivity  gradient  [Atkinson  and 
Hutchison,  1978].  The  polar  cap  potential  drop  is  about  60 
kV,  in  good  agreement  with  empirical  values  for  the  prevail¬ 
ing  A'p-index  [e.g.,  Foster  et  al.,  1986;  Heppner  and  May¬ 
nard,  1987],  and  directed  roughly  between  0200-0300  MLT 
and  1500  MLT. 

Figure  7  (bottom)  shows  the  model  and  calculated  poten¬ 
tial  profile  along  the  ionospheric  projection  of  the  Viking 
orbit.  This  serves  to  demonstrate  that  the  model  results  are 
not  only  qualitatively  but  also  quantitatively  in  good  agree¬ 
ment  with  the  observations.  Deviations  occur,  however,  in 
two  time  intervals,  most  pronounced  around  0830  DT  and 
less  clear  around  0720  UT,  where  the  high-altitude  poten¬ 
tial  (Viking)  has  minima  which  are  smaller  or  absent  in  the 
ionospheric  (model)  potential.  It  is  evident  from  the  mea¬ 
sured  field-aligned  current  profile  (top)  that  intense  upward 
currents  are  present  in  these  regions.  Furthermore,  these 
regions  show  evidence  of  upward  acceleration  of  ions  as  in¬ 
dicated  by  the  inserted  shaded  areas.  This  is  exemplified 
in  Figure  8,  which  shows  pronounced  upward  ion  beams  in 
the  two  highest-energy  channels  (0.52  keV  and  1.25  keV) 
between  0720  and  0723  LIT.  In  the  downward  current  re¬ 
gion  between  0718  UT  and  0720  UT,  weak  ion  conics  can 
be  seen  at  the  low-energy  channel  of  60  eV.  The  electron 
data  (not  shown  here)  show  only  weak  acceleration  signa¬ 
tures  end  are  thus  consistent  with  the  upward  ion  beams 
observed  essentially  above  the  acceleration  region.  Another 
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Fig  7.  Quantitative  point-by-point  comparison  between  model  and  observed  potential  profiles  (middle)  and 
between  model  and  observed  field-aligned  currents  (lop)  along  the  ionospheric  projection  of  the  Viking  orbit.  The 
shaded  areas  in  the  bottom  indicate  regions  of  upward  accelerated  ions  and  intense  AKR 
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Fig  R  Timr  s<Tirs  of  ion  fluxfs  at  four  sflrcted  rnfrgifs  (0  06.  0  22,  0  52,  1.25  k-V)  obtained  by  the  ion  detectors 
(PISPI/2J  on  Viking  The  ion  fluxes  are  shown  on  a  linear  scale  in  relative  units.  The  basic  flux  level  in  units  of 
(keV  cm^  5  sr)“*  is  given  at  the  right  side  of  each  panel  The  pitch  angle  is  Indicated  in  the  lowest  panel  The 
ion  fluxes  are  seen  to  be  concentrated  around  180°  pitch  angle  (upward). 


piece  of  evidence  that  Viking  passed  above  an  auroral  ac¬ 
celeration  region  centered  around  0720  UT  is  given  by  the 
high- irequency  wave  spectrogram  in  Figure  9.  The  elec¬ 
tric  component  is  shown  here  for  the  time  interval  between 


0717:20  UT  and  0722:10  UT.  Intense  auroral  kilometric  ra¬ 
diation  (AKR)  emissions  centered  around  0719:40  UT  with 
low-frequency  cutoffs  between  350  and  400  kHz  are  seen 
at  the  top  well  above  the  electron  gyrofrequency  indicated 


Viking  High-frequency  Electric  Field 


AKR 


7:19 


9  Electric  field  observed  by  the  high-frequency  wuve  experiment  on  Viking  during  the  afternoon  auroral 
oval  crossing  around  0720  UT.  Intense  AKR  emissions  are  seen  at  the  top  well  above  the  electron  gyrofrequency 
indicated  by  the  solid  line  around  100  kH*  Impulsive  hiss  emissions  are  seen  at  the  bottom,  shortly  before  0720 
UT 
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midnight.  These  currents  can  be  compared  with  currents 
that  could  be  inferred  from  simultaneous  ground-based  mag¬ 
netometer  observations.  The  Joule  power  dissipation  is  seen 
to  be  concentrated  into  two  regions:  (1)  the  midnight  sector 
and  in  particular  the  bright  auroral  structures  there,  and  (2) 
the  prenoon  sector,  characterized  by  relatively  intense  auro¬ 
ral  structures  seen  by  the  UV  imager  and  relatively  intense 
field-aligned  currents  as  measured  by  both  the  DMSP  F7 
and  Viking. 

A  simplified  picture  of  the  magnetospheric  convection  is 
illustrated  in  Figure  11  (right),  which  is  a  projection  of  the 
deduced  ionospheric  two-cell  convection  pattern  (left)  to  the 
equatorial  plane  using  a  dipole  magnetic  field.  The  location 
of  the  last  closed  equipotential  contour  (heavy  line)  is  seen  to 
correspond  relatively  well  with  that  of  the  average  plasma- 
pause  boundary  (dotted  line)  as  given  by  Carpenter  [1966]. 
(Since  the  shape  of  the  plasmapause  varies,  a  detailed  cor¬ 
respondence  is  not  to  be  expected.) 


Fig  10  The  ionospheric  current  vectors,  plotted  every  2°  of 
invariant  latitude  and  every  hour  of  magnetic  local  time  (top), 
and  Joule  power  dissipation  (bottom)  corresponding  to  the  poten¬ 
tial  distribution  shown  in  Figure  6  and  conductivity  distribution 
shown  in  Figure  2  (top). 


by  the  solid  line  at  around  100  kHz.  Note  also  the  impul¬ 
sive  character  of  the  low-frequency  hiss  emissions  around 
0719:55  UT,  which  is  characteristic  for  the  crossing  of  an 
arc  structure. 

To  summarize,  the  model  output  potential  distribution 
agrees  both  qualitatively  and  quantitatively  with  the  Viking 
observations,  except  in  two  regions  where  expected  devi¬ 
ations  occur,  qualitatively  consistent  with  upward  parallel 
electric  fields. 

As  a  further  illustration  of  the  wide  range  of  results  that 
can  be  obtained  with  this  model,  Figure  10  shows  the  iono¬ 
spheric  current  (top)  and  the  Joule  power  dissipation  (bot¬ 
tom)  associated  with  this  event  and  calculated  from  the 
model  electric  field  and  conductivity  data  using  equations 
(1)  and  (11),  respectively.  Note  the  intense  westward  elec¬ 
trojet  associated  with  the  auroral  activity  close  to  magnetic 


6.  Discussion 

UV  imaging  of  the  aurora  from  polar-orbiting  satellites 
such  as  Dynamics  Explorer  1  and  Viking  represents  an  im¬ 
portant  new  input  to  the  study  of  auroral  electrodynamics 
on  a  global  scale.  The  images  of  the  auroral  oval  provide 
not  only  a  natural  reference  frame  for  the  other  observa¬ 
tions  but  also  valuable  information  on  various  electrody¬ 
namical  parameters.  It  is,  for  example,  evident  that  the 
observed  variations  in  the  UV  emission  intensities  also  re¬ 
flect  variations  in  the  ionization  and  thus  the  ionospheric 
conductivity.  This  dependence  has  been  used  by,  for  exam¬ 
ple,  Craven  et  ai.  [1983]  and  KamJde  et  aJ.  [1986]  to  infer 
“instantaneous”  conductivity  distributions  from  DE  1  au¬ 
roral  images.  A  quantitative  estimate  involves,  however,  a 
number  of  uncertainties,  as  has  been  discussed  in  section  2. 

It  is  also  clear  that  active  regions  such  as  bright  auroral 
structures  are  associated  not  only  with  enhanced  conductiv¬ 
ities  but  also  with  enhanced  upward  field-aligned  currents. 
It  b,  however,  very  unclear  how  the  information  of  the  fine- 
scale  structure  seen  in  the  weak  and  diffuse  background  oval 
could  be  used  to  better  reconstruct  the  associated  field- 
aligned  ciuient  dbtribution.  Therefore,  it  makes  no  sense 
to  represent  these  fine-scale  variations  in  the  model  input 
data  unless  it  could  be  done  not  only  for  the  conductivity 
but  abo  for  the  field-aligned  current. 

For  these  reasons  we  have  chosen  an  alternative  and  qual¬ 
itative  way  of  utilizing  the  UV  image  information,  namely, 
to  obtain  a  reference  frame  for  the  modeling  and  to  identify 
the  locations  and  extent  of  the  major  intense  auroral  struc¬ 
tures,  which  are  modeled  by  enhancements  in  the  preexisting 
quiet  time  field-aligned  current  [lijinta  and  Potemra,  1976a, 
1978]  and  conductivity  dbtributions,  here  representative  of 
the  background  auroral  oval. 

An  important  and  new  feature  of  the  model  presented 
here  b  that  the  field-aligned  current  and  conductivity  pat¬ 
terns  are  consbtent  with  each  other  in  the  sense  that  the 
conductivity  peaks  in  regions  of  upward  field-aligned  cur¬ 
rents  associated  with  the  discrete  aurora.  To  determine  the 
exact  form  of  such  a  correspondence  b  a  very  complicated 
problem  since  the  precipitating  particles  creating  enhanced 
ionisation  are  typically  more  energetic  than  those  carrying 
the  major  part  of  the  upward  field-aligned  currents.  Qual¬ 
itatively,  such  a  correspondence  b,  however,  well  confirmed 


I4.49() 


MaRKLUND  BT  AL.:  SNAPSHOTS  OP  HlOH-LATITUtiB  ELECTRODYNAMICS 


0 


Inertial 


IONOSPHERE  MAGNETOSPHERE 

Equatorial  plane 


Fig.  11.  The  deduced  two-cell  convection  pattern  (left)  and  the  corresponding  magnetospheric  convection  (right) 
resulting  from  a  sunple  projection  to  the  equatorial  plane  using  a  dipole  magnetic  field.  For  comparison  the  average 
plasmapause  from  Carpenter  [1966]  has  been  inserted  (dotted  line). 


from  observations  over  auroral  arcs.  Lyons  et  ai.  [1979] 
found  that  when  a  field-aligned  acceleration  voltage  (V) 
could  be  inferred  from  the  electron  spectra,  the  net  down¬ 
ward  electron  energy  flux  generally  varied  as  V’  while  the 
net  upward  field-aligned  current  was  generally  proportional 
to  V.  Since  Sp  is  roughly  proportional  to  the  square  root 
of  the  electron  energy  flux  (see  section  2  above),  the  rela¬ 
tions  above  imply  a  direct  proportionality  between  Sp  and 
yi|.  A  correspondence  between  these  two  parameters  is  also 
evident  from  Figures  3  and  5,  which  justifies  our  choice  of 
the  linear  relationship  described  by  equation  (10). 

The  new  method  described  in  this  paper  is  still  under 
development,  and  refinements  of  some  of  the  details  in  the 
model  are  to  be  expected  as  soon  as  we  have  gained  more 
experience  from  several  different  events.  The  criteria  used 
in  the  selection  of  a  suitable  event  for  a  first  study  were 
that  the  entire  auroral  oval  should  be  visible  to  the  imager; 
that  at  least  one  large-scale  auroral  structure  should  exist 
ill  some  local  time  sector  of  the  oval;  and  that  the  substorm 
activity  should  be  relatively  stable  during  the  time  period 
for  which  the  UV  image  is  representative.  The  event  chosen 
here  from  Viking  orbit  360  was  found  to  meet  these  various 
requirements  reasonably  well.  The  requirement  of  intense 
aurora  implies,  by  definition,  also  substorm  activity,  which, 
as  demonstrated  by  the  AE  index,  was  high  but  relatively 
stable  (500  nT)  during  the  time  interval  of  27  min  for  which 
the  UV  image  is  assumed  to  be  representative. 

For  the  daytide  auroral  region  we  do  not  see  any  problem 
with  this  assumption,  since  all  information  used  to  calibrate 
the  model  input  data  (from  the  Viking  afternoon  auroral 
oval  crossing  and  the  DMSP  F7  prenoon  auroral  oval  cross¬ 
ing)  was  obtained  within  15  min.  Since  the  triple  field- 
aligned  current  sheets  iii  the  prenoon  sector,  as  observed 


by  DMSP  F7  around  0739  UT  and  Viking  about  1  hour 
later,  had  very  similar  characteristics  (in  terms  of  magni¬ 
tudes,  widths,  etc.),  the  dayside  activity  might  have  been 
relatively  stable  even  during  the  entire  auroral  oval  transit 
time  by  Viking.  This  is  also  supported  by  the  good  agree¬ 
ment  obtained  between  the  model  potential  distribution  and 
that  inferred  from  the  measurements  (cf.  Figures  6  and  7). 

For  the  nightside  auroral  oval  the  requirement  of  roughly 
unchanged  auroral  conditions  relative  to  those  of  the  UV 
image  during  27  min  cannot  be  verified  from  the  data  set 
available  here.  We  have  assumed  this  to  be  the  case,  sup¬ 
ported  by  the  observations  of  intense  particle  fluxes  and 
field-aligned  currents  in  the  region  where  the  bright  auroral 
intensification  was  observed  by  the  UV  imager  half  an  hour 
earlier. 

The  global  field-aligned  current  and  conductivity  distri¬ 
butions  presented  in  Figure  2  and  the  corresponding  equipo- 
tential  pattern  presented  in  Figure  6  are  believed  to  rep¬ 
resent  the  prevailing  electrodynamical  situation  associated 
with  the  auroral  situation  illustrated  by  Figure  1  .  The  qual¬ 
itative  agreement  between  the  model  convection  pattern  smd 
the  E  X  B  drift  vectors  and  the  good  quantitative  agreement 
between  the  model  and  calculated  potential  profile  along  the 
Viking  orbit  are  a  consequence  of  the  careful  selection  of  re¬ 
alistic  global  input  data  being  calibrated  to  the  Viking  and 
DMSP  F7  data.  The  asymmetry  between  the  intense  and 
crescent -shaped  dawn  convection  cell  (associated  with  the 
intense  auroral  activity  and  field-aligned  currents)  and  the 
much  weaker  dusk  convection  ceU  shown  in  Figure  6  is  con¬ 
sistent  with  earlier  observations  and  theoretical  predictions 
for  positive  IMF  By  conditions  [Heppner,  1972;  Moser  et  al., 
1974;  Crooker,  1979). 

It  is  interesting  that  the  “instantaneous”  potential  pattern 
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calculated  here  from  model  input  conductivity  and  field- 
aligned  current  distributions  is  quite  similar  to  the  aver¬ 
age  BC  pattern  of  Heppner  and  Maynard  [1987]  based  on 
a  large  statistical  material  of  DE  2  electric  field  observa¬ 
tions  and  representative  for  positive  IMF  By  conditions  and 
3+  <  Kp  <  4  (representative  also  for  Viking  orbit  360). 
The  average  direction  of  the  polar  cap  electric  field  between 
0300-0400  MLT  and  1500-1600  MLT,  the  polar  cap  poten¬ 
tial  drop  of  60  kV,  the  dominant  morning  convection  ceil, 
and  intrusion  of  the  morning  cell  potential  contours  toward 
west  in  the  nightside  auroral  oval  and  toward  east  in  the 
dayside  auroral  oval  are  features  which  are  characteristic  for 
both  the  "instantaneous''  and  the  statistical  pattern. 

In  order  to  combine  data  from  various  satellites  follow¬ 
ing  the  procedure  outlined  here  a  constructive  assumption 
is  that  there  is  a  perfect  coupling  between  the  electrody¬ 
namical  parameters  in  the  ionosphere  and  at  satellite  alti¬ 
tude.  The  good  agreement  between  the  measurements  made 
at  the  two  satellite  altitudes,  as  well  as  between  the  mea¬ 
surements  and  the  model  input  and  output  data,  generally 
supports  this  assumption.  As  a  result  of  using  the  assump¬ 
tion  E  •  B  =  0  we  found,  however,  in  the  regions  of  in¬ 
tense  upward  field-aligned  currents,  discrepancies  between 
the  model  and  calculated  potentials  (cf.  Figures  6  and  7) 
indicative  of  upward  parallel  electric  fields.  U'ei'mer  et  af. 
[1985]  compared  the  electric  field  signatures  at  high  and  low 
altitudes  as  measured  by  the  Dynamics  Explorer  1  and  2 
satellites.  In  agreement  with  the  findings  here  they  found 
that  the  large-scale  characteristics  of  the  electric  fields  at 
the  two  altitudes  were  very  similar  but  that  the  small-scale 
electric  field  in  the  auroral  zone  had  larger  magnitudes  at 
higher  altitudes  consistent  with  the  existence  of  parallel  elec¬ 
tric  fields. 

The  presence  of  upward  accelerated  ions  and  intense  AKR 
emissions  (0720  UT)  well  above  the  local  electron  gytofre- 
quency  and  absence  of,  or  only  weak,  signatures  of  electron 
acceleration  in  these  regions,  all  indicate  that  Viking  passed 
above  or  crossed  the  upper  part  of  the  acceleration  regions. 
For  the  two  regions  of  upward  field-aligned  currents  indi¬ 
cated  in  Figure  7  the  potential  drop  amounts  to  about  5-10 
kV,  which  is  within  the  typical  range  of  acceleration  energies 
for  the  auroral  particles. 

The  technique  presented  here  is  in  principle  applicable 
also  to  more  disturbed  and  thus  time-dependent  events. 
Several  UV  images  are  then  needed  to  follow  the  tempo¬ 
ral  evolution  of  the  aurora  as  the  satellite  crosses  over  it.  In 
this  kind  of  time-dependent  modeling,  calibration  to  simul¬ 
taneous  measurements  by  several  spacecraft  is  particularly 
useful. 

7.  Summary 

Simultaneous  observations  by  the  Viking  and  DMS?  F7 
satellites  have  been  used  in  a  new  method  to  obtain  snap¬ 
shot  pictures  of  the  high-latitude  electrodynamics  associ¬ 
ated  with  a  specific  auroral  situation  viewed  by  the  UV 
imager  experiment  on  Viking  during  orbit  360.  The  UV 
image  was  used  to  locate  the  regions  of  active  intense  auro¬ 
ras,  which  were  represented  by  enhanced  conductivities  and 
upward  field-aligned  currents  in  the  model  input  data.  The 
global  input  distributions  were  calibrated  so  as  to  be  quan¬ 
titatively  consistent  with  the  conductivity  and  field-aligned 
current  profiles  calculated  from  the  particle  and  magnetic 
field  data  along  the  two  satellite  orbits.  These  distribu¬ 


tions  were  used  together  with  Ohm’s  law  to  solve  for  the 
ionospheric  potential  distribution,  which  was  finally  tested 
against  the  electric  field  observed  by  Viking.  The  results 
can  be  summarized  as  follows; 

1.  The  model  input  data  (E  and  y||)  are  qualitatively 
consistent  with  the  auroral  imager  data  and  quantitatively 
consistent  with  the  Viking  and  DMSP  F7  observations.  An 
important  feature  in  the  model  is  a  current  (yji  )-conductivity 
coupling  term. 

2.  There  is  a  good  agreement  both  qualitatively  between 
the  output  potential  distribution  and  the  observed  E  x  B 
drift  and  between  the  two  corresponding  potential  profiles 
along  the  Viking  orbit.  The  “instantaneous”  potential  dis¬ 
tribution  is,  further,  quite  similar  to  the  empirical  BC  model 
of  Heppner  and  Maynard  [1987]  based  on  DE  2  electric  field 
data  ftom  periods  of  positive  IMF  By  and  3+  <  Kp  <  4-, 
which  is  the  same  conditions  as  for  the  event  studied  here. 

3.  Potential  differences  between  the  Viking  altitude  and 
the  ionosphere  indicative  of  upward  parallel  electric  fields  are 
found  for  two  time  periods  corresponding  to  Viking  crossings 
above  acceleration  regions  further  characterized  by  upgo> 
ing  ions,  intense  upward  field-aligned  currents,  and  intense 
AKR  emissions. 

4.  To  illustrate  the  variety  of  model  results  that  can 
be  obtained,  the  ionospheric  currents,  Joule  power  dissipa¬ 
tion,  and  magnetospheric  convection  representative  for  this 
event  have  been  estimated  using  the  model  conductivity  and 
potential  distributions.  An  equatorial  plane  projection  of 
the  potential  pattern  shows  for  the  last  closed  equipotential 
contour  a  good  agreement  with  the  average  location  of  the 
piasmapause. 
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